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1 . INTRODUCTION 


1.1 General 

This document presents the final report on the Shuttle/TDRSS 
Communications Systei.. Performance Analysis performed for NASA Johnson 
Space Center under Contract NAS 9-15799 directed by William Teasdale 
and Sid Novosad. It represents a portion of the wn»k accomplished 
during the period April 20, 1979 through January 24, 1980. 

The general objectives of the overall contract are the following: 

(1) T o modify and refine the existing Shuttle/TDRSS link simulation 
programs to model the post-RFI TDRS hardware and to provide the 
ability to evaluate the link performance degradation due to 

RFI effects. This analysis had to include the S-band synchroniza- 
tion vulnerability to the TDRSS RFI environment. 

(2) To use the refined link models to determine, evaluate and assess 
expected link performance that will result from not meeting each 
TDRSS user constraint; individually, and in various combinations. 

The effects and system performance implications of the latest user 
constraint values had to be evaluated. The refined results of 

this task were to be used by JSC as a technical base for negotiating 
interface parameters between Shuttle and TDRSS. 

(3) To continue the review of ESTL Task 501 TDRSS system configuration 
(using flight equivalent TDRSS hardware) and recommend changes 
an/or modifications as required. LinCom continued its support 

of ESTL Task 501 TDRSS system performance prediction capability 
by providing math modeling and computer models necessary to provide 
valid system level performance predictions for forward and return 
link certification testing. 
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(4) To identify test requirements to evaluate user-constraints for 
the refined TORSS model and update Shuttle Ku-band and S-band 
system characteristics. 

In what follows an overall description of the contractual effort 
and a brief summary of the results is given. This is followed by backup 
material which includes simulation data and analyses from which our 
summary results and recommendations have been drawn. 

1.1.1 Final Report Contents 

This report addresses and documents LinCom's findings on the 
task statements detailed in the Statement of Work. All task statements 
pertain to the performance prediction of Shuttle forward and return 
links through the Tracking and Data Relay Satellite System (TDRSS). 

Chapter 2 presents user constraint sensitivity data for Shuttle/ 
TDRSS S-band and Ku-band return links. These simulation results are 
based on refined TDRS model and up-to-date user constraint values. 

They apply to links without RFI problems. The current system 
parameter values are also listed. 

Chapter 3 contains a comprehensive performance analysis for the 
Shuttle/TDRSS ground station carrier and timing recovery circuits. 

It contains performance data for the current system parameters and 
link budget for both Ku-band links and S-band links. The effects of 
RFI are not considered in this analysis. 

Chapter 4 discusses the effects of RFI on link performance. The 
RFI environment is described in detail and the RFI-related hardware 
changes in the TDRS and TDRSS ground station are indicated. The 
analytical models used to characterize the bit error rate and 
synchronization performance are described. Preliminary results for 
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the BER and synchronization vulnerability to RFI are also included. 

In Chapter 5 the problem of simulating the RFI effects with 
ESTL's flight equivalent TDRSS hardware is addressed. The RFI pulse 
characteristics which might have to be considered in the design of 
an RFI test generator are discussed and a simple implementation based 
on unclassified RFI characteristics is presented. 

Chapter 6 comprises a number of small studies performed under 
this contract in response to specific requests by JSC. 

1 .2 Summary and Recommendations 

1.2.1 Shuttle S-band Return Link Performance Study 

The performance degradation of the Shuttle/TDRSS return link 

due to Shuttle signal imperfections has been evaluated for the 

expected Shuttle distortion values based on current link budgets 

and TDRS/grour.d station hardware data. The results show a loss 

relative to ideal BPSK of 1.65 dB at the nominal bit error 
-4 

probability of 10 . The TDRS and ground station contribute 
.9 dB to this loss and the Shuttle the remaining .75 dB. This is 
illustrated in Fig. 1.1. 

The single biggest contributor to the Shuttle degradation 
is the phase noise which is specified as 10° rms in the frequency 
interval 0 to 270 Hz. We recommend to separate the phase noise 
into components lying inside and outside the carrier tracking loop 
bandwidth, respectively, since these components affect the bit 
error rate performance differently. Such a refinement of the 
specification could greatly improve the accuracy of the predicted 
performance results. Since our present prediction is based on 
worst-case assumptions the overall degradation would be reduced. 
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Fig. 1.1. Performance of Shuttle S-Band Pveturn Link, 
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The effects of RF I both on the bit error rate and the synchronization 
performance of the Shuttle/TDRSS S-band return link was evaluated. 

For the TDRS West (which sees a less severe RFl environment than 
the TDRS East) and the S-band frequency, 2217.5 MHz the TER was 
found to be degraded by .2 dB for 4 degrees off-pointing, by .6 dB 
for 1.5 degrees off-pointing and by almost 6 dB when the Shuttle is 
directly over the RFI zone. This is illustrated in Fig. 1.2. 



- 5 - 





- 6 - 




JlnCc 


om 


Table 1-2. Timing Recovery Performance. 


Channel 

Data Rate 

V N 0 

RMS Jitter (Symbol Fraction) 

Ku-Band 
Channel 1 

192 kb/s 

10 dB 

.056% 



f 


K-Band 
Channel 2 

600 kb/s 

11 dB 

.15% 

Ku-Band 
Channel 2 

10 Mb/s 

4 dB 

2.2% 

S-Band 

96 kb/s 

4 dB 

.68% 
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1.2.5 ESTL Hardware RFI Testinc 


The question of designing an RFI test generator for use in conju- 
tion with the ESTL TDRSS flight equivalent Hardware is addressed and 
a simple test generator implementation is proposed. This design is 
based on unclassified data bases describing the RFI statistics. It 
is recommended to review the more detailed classified data base prior 
to the final hardware design to insure that all pertinent RFI pulse 
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2. LinCsim PREDICTIONS OF SHUTTLE/TDRSS RETURN LINK BER PERFORMANCE 

2.1 Introduction 

The Shuttle return link communications hardware does not meet all 
the specifications set forth in the User Constraint Table of the TDRSS 
Users' Guide [1], Since the TDRSS performance predictions are based 
on these interface parameter values they may not be applicable to the 
Shuttle links. In this chapter the performance sensitivity to these 
user constraint values is documented. The design point is based on the 
latest hardware data and user constraint values available. All data 
apply for a non-RFI environment. 

Section 2.2 presents sensitivity data for the Shuttle/TDRSS S-band 
return link. Mode 2. Since this is the high rate mode it is more 
susceptible to most of the distortion effects and the resulting sensitivities 
are therefore an upper bound on the degradations in the low-rate mode. 

Section 2.3 presents the results for all three channels of the Ku-band 
return link. Mode 1. 

2.2 Performance Prediction for Shuttle S-band Return Link 

This section presents the sensitivity of the Shuttle S-band return 
link BER performance to variations of tne user constraint values 
based on current link budgets, nominal user constraint value estimates 
and TDRS/grour.d station hardware data. 

The link characteristics assumed are summarized in Table 2.1. 

The nominal values of the user constraints are listed in Table 2.2. 

They agree with the expected Shuttle values obtained from JSC. The 
link budget is reproduced in Table 2.3. It is based on the budgets 
contained in Ref. 2 with some updates obtained from Dr Kwei Tu. 

Table 2.4 lists the user constraints whose effect on performance wjs 
studied. 
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Table 2.1. Shuttle S-band Return Link Characteristics 
Used for LinCsim 


Data Rate 
Data Format 
Coding 

Carrier 

Design Error Rate 


192 Kbps 

BPSK, eiphase. Unspread 

Rate 1/3 Convolutional 
Code 

Noncoherent with Forward 
Link 
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Table 2.2. Nominal User Constraint Values Used for 
LinCsim 


Data Bit Jitter (3o) 


Modulator Phase Imbalance 

Modulator Gain Imbalance 

Data Asymmetry 

Phase Nonlinearity 

Gain Flatness 

Gain Slope 

AM/PM 

3 dB Bandwidth 


Phase Noise 


1 Hz - 10 Hz 
10 Hz - 100 Hz 
100 Hz - 1 kHz 
1 kHz - 6 MHz 


.4 dB peak 
.1 dB/MHz 
14°/dB 
100 MHz 
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Table 2. 3. Shuttle S-Band Return Link Power Budget 


Shuttle-to-TDRS Link 
Shuttle EiRP 
Space Loss 
Polarization Loss 
TORS G/T 

TDRS-to-Ground Link 

TDRS TWT max Output Power 
TWT Output Backoff 
SSA Power Allocation 
TORS Hardware Losses 
TDRS Antenna Gain 
Pointing Loss 
Space Loss 
Polarization Loss 


16.7 dBW 
192.1 dB 
.5 dB 

9.55 dB/°K 

13.4 dBW 
2 dB 
<10 dB 
2.46 .iB 
45.9 dB 
.65 dB 

207.7 dB 
.1 dB 


Atmospheric Loss 
Ground Station G/T 


.8 dB 

40.3 dB/°K 
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Table 2.4. Parameters ' tudied. 


Modulator Gain Imbalance 

Modulator Phase Imbalance 

Data Asymmetry 

Data Bit Jitter 

Data Static Timing Offset 

XTR AM/AM 

XTR AM/PM 

Static Phase Error 

XTR Gain Slope 

XTR Gain Flatness 

XTR Phase Nonlinearity 

Phase Noise R.M.S. 
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The results are shown in two different forms. The BER curves 
show the bit error probability as a function of Shuttle-to-TDRS link 
carrier-to-noise ratio (CNR) variation around the nominal link budget 
of Table 2.3. The horizontal distance between the bit error rate curve 
and the design point shown represents the margin in carrier-to-noise 
ratio which can be allocated to the various subsystems for hardware 
degradations. The sensitivity curves show the increase in the Shuttle- 
to-TDRSS CNR needed to offset the performance degradation (relative to 
the nominal performance) due to the variation of a single parameter at the 

j* 

design error rate of 10 . 

The error rate curve in Fig. 2.1 represents the BER performance of 
a Shuttle transponder transmitting a perfect signal (i.e., all 
user constraints are set to zero). The CNR loss shown (.9 dB) 

can be attributed to the TDRS and the ground station. This BER curve 
is reproduced on all other BER plots and labeled "Perfect User." The 
horizontal distance between this curve and one of the other BER curves 
represents the CNR loss due to the combined effect of all the user 
constraints. For the nominal conditions this loss amounts to .75 dB. 

Fig. 2.20 shows that by far the biggest contribution of degradation 
comes from the phase noise which is specified as 10° rms in the 0 to 
270 Hz frequency range. As a worst-case assumption this noise power 
was concentrated in the 100-270 Hz range (i.e. outside the tracking 
loop bandwidth). A more accurate performance prediction could be 
made if the noise power inside and outside the tracking loop bandwidth 
were known separately. 


Tables 2.5 and 2.6 contain a complete list of all nominal system 
parameters used in LinCsim, including those pertaining to the TDRS and 
TDRSS ground station. 
















































































Table 2.6. Complete List of Simulation Parameters 


PARAMETER 


PRESENT 

VALUE 


WHERE 

OBTAINED 


DATE 

OBTAINED 


OLD VALUE 


SOURCE 
OLD VAIUE 


MR EIRP: UNCOOEO 
DATA 


AIR EIRP: CODEO 
DATA 


SPACE LOSS 


POLARIZATION LOSS 


ATMOSPHERIC LOSS 


TORS G/T 


TORS TWT MAX 
OUTPUT POWER 


TORS TWI OUTPUT 
BACKOFF 


POWER ALLOCATION 


CHANNEL POWER 


TORS HARDWARE LOSS 


TDRS ANTENNA GAIN 


POINTING LOSS 


SPACE LOSS 


POLARIZATION LOSS 


ATMOSPHERIC LOSS 


RX ANIL SNA G/T 


XT OSC .COHERENT 
TURNAROUND 


|'j.S5(po"i)na 


13.4 dB 


1.4 dB 


2.45 d5 


45.99 dB 


0-65 tjg 


HZ to 1 kHz 


z 


XT OSC NONCOHERENT 
TURNAROUND 


to 10 Hz 


Z toTiHZ 


Z 


TDRS OSCILLATORS 


RX OSCILLATORS 


RX CARRIER TRACKING 
LOOP 


STATIC PHASE 


Masai 


PHASE NOISE BUDGET 


T9. 2 dBW 


8 

3 

dB/'K 

12 

5 

dBW 

2 

0 

dBW 

m 

y 

oBW 

0 

6 

dBW 

1 

7 

dB 

41 

3 

dB 

0 

7 

dB 

207 

7 

dB 

0 

1 

dB 

0 

8 

dB 

40 

3 

dB 


ISS Circuit Margin,V77 




ISS Circuit M.?rgin,V77 


U Proposal ,7/76 


U Froposal .7/76 


!U Proposal , 7/76 


WU Proposal , 7/76 


WU Proposal, 7/ 7b 


Proposal , 7/76 


WU Proposal, 7/76 


III Proposal , 7/76 


J Proposal, 7/76 


S-B05-T 


S-80S-1 


inCom Estimate 


HTOTSKH 


squaring LOSS 


LinCom Estimate 
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Figure 2 . 2 . BE R Plot for Gain Imbalance. 
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Figure 2.3. Sensitivity Curve. 
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Figure ?,6.BER Plot for Data Asyrunetry. 
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Figure 2 _8. BER Plot for Bit Jitter. 
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Figure 2.9 .Sensitivity Curve. 
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Figure 2.1.2. BEF. Plot for AM/PM 
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Figure 2.13. BER Plot for Static Phase Error. 
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Figure 2.1.4. Sensitivity Curves. 
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Figure 2.15. BER Plot for Gain Slope. 
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Figure 2.16. BER Plot for Gain Flatness. 


71 * 00 *.; 3 








Figure 2 .18. Sensitivity Curves . 
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Figure 2 . 21 . BER Plot for High Frequency Noise. 
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Figure 2. 22. Sensitivity Curves for Phase Noise. 
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2.3 Performance Prediction for Snuttle Ku-Band Return Link, Mode 1 


This section presents the sensitivity of the Shuttle/TDRSS 
Ku-band return link BER performance to variations of the user 
constraint values. The all-digital mode 1 is considered and results 
are given for all three data channels. The link characteristics 
assumed are summarized in Table 2.7. The nominal values of the user 
constraints are listed in Table 2.8. They agree with the expected 
Shuttle values whereever such data were available from JSC. The 
link dubget is reproduced in Table 2.9; it is based on Ref. 2. 

Table 2.10 lists the user constraints whose effect on performance was 
studied. 

The results are shown in two different forms. The BER plots show 
the bit error probability as a function of the Shuttle- to-TDRS link 
carrier-to-noise ratio (CNR) variation around the nominal link budget 
of Table 2.9. The horizontal distance between the bit error rate curve 
and the design point shown represents the margin in carrier-to-noise 
ratio which can be allocated to the various subsystems for hardware 
degradations. The sensitivity curves show the increase in Shuttle-to- 
TDRSS CNR needed to offset the performance degradation (relative to 
the nominal performance) due to the variation of a single parameter at 
the design error rate of 10~®. 

The error rate curves in Figs. 2 .23, 2 .40 and 2 .57 represent the 
BER performance of channels 1, 2, and 3, respectively, for a Shuttle 
transponder transmitting a perfect signal (i.e. all user constraints 
are set to zero). The CNR loss shown (Table 2. 11) can be attributed to 
the TORS and the ground station. This BER curve is reproduced on all 
other BER plots and labeled "Perfect Shuttle." The horizontal distance 
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between this curve and one of the other 8ER curves represents the 
CNR loss due to the combined effect of all the user constraints. 

For the nominal conditions this loss is shown in Table 4.11. Table 
4.12 lists the major contributors to this loss for each of the three 
channels. 

Table 4.13 contains a complete list of all nominal system 
parameters used in LinCsim, including those pertaining to the TDRS 
and TDRSS ground station. 
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Table 2.7 . Shuttle Ku-Band Return Link Characteristics 
Used for LinCsim. 


Data Rate 

Channel 

1 

192 kb/sec 


Channel 

2 

2 Mb/sec 


Channel 

3 

50 Mb/sec 

Date. Format 

Channel 

1 

Bi -phase 


Channel 

2 

NRZ 


Channel 

3 

NRZ 

Coding 

Channel 

1 

none 


Channel 

2 

none 


Channel 

3 

rate 1/2 


Carrier Noncoherent with Forward Link 

Square Wave 


Subcarrier 
Design Erro' Rate 
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Table2.8. Nominal User Constraint Vafues Used for LinCsim 


Relative Phase Between 
I and Q Channels 2° 

Data Asymmetry 10% 

Modulator Gain Imbalance .5 dB 

Gain Flatness .3 dB 

Gain Slope .01 dB/MHz 

Phase Nonlinearity 7° 

AM/PM 5°/dB 

Data Bit Jitter .01% rms 

3 dB Bandwidth 200 MHz 

Carrier Phase Noise 

100 Hz - 1 kHz 17° rms 

1 k Hz - 150 MHz 3° rms 

Subcarrier Phase Noise 1° rms 
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Table 2. 9. Shuttle Ku-Band Return Link Power Budget. 


Shuttle-to-TDRS 


Shuttle EIRP 

48.1 dBW 

Space Loss 

208.5 dB 

Polarization Loss 

.3 dB 

TDRS G/T 

22.6 dB/°K 

i-to-Ground Link 

TDRS TWT max Output Power 

12.5 dBW 

TWT Output Backoff 

.5 dB 

TDRS Hardware Losses 

1.2 dB 

TDRS Antenna Gain 

41.3 dB 

Pointing Loss 

.7 dB 

Space Loss 

207.7 dB 

Polarization Loss 

.1 dB 

Atmospheric Loss 

.8 dB 

Ground Station G/T 

40.3 dB/°K 
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Table 2.10. Parameters Studied. 

Carrier Power Split 

Modulator Gain Imbalance 

Relative Phase Between I and Q Channels 

Data Asymmetry 

Data Bit Jitter 

Transmitter AM/AM 

Transmitter AM/PM 

Carrier Phase Noise 

Gain Flatness 

Gain Slope 

Phase Nonlinearity 
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Table 2.12.Major Contributors to CNR Loss. 


Channel 1 : 

Modulator Gain Imbalance 
AM/ AM 

Data Asymmetry 

Subcarrier Power Split Error 


Channel 2 : 

Modulator Gain Imbalance 
Carrier I/Q Relative Phase 
AM/AM 
AM/PM 

Data Asymmetry 


Channel 3: 

Modulator Gain Imbalance 
Carrier I/Q Relative Phase 
Data Asymmetry 
AM/PM 

Phase Nonlinearity 
Gain Flatness 
Gain Slope 
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Table 2 .13 (cont'd) 


PARAMETER 

PRESENT 

VALUE 

WHERE 

OBTAINED 

DATE 

OBTAINEO 

OLO VALUE 

SOURCt 
OLD VALUE 


1 LINK BuoC'ET 1 1 

* 




‘ 


XTR EIRP 




43.1 dBW 

SSRF Circuit Margin, 
1/77 

SPACE LOSS 




208.5 dB 

SSRF Circuit Margin, 
1/77 . 

POLARIZATION LOSS 




-6.3 dB 

Circuit Margin, 

ATMOSPHERIC LOSS 




0 dB 

SSRF CircuitMargin 1/77 

TORS G/T 

* 



22.6 dB/°K 

SSRF Circuit Marqin T /7' 

TORS TUT MAX 
OUTPUT POWER 


- 


12.5 dB 

WU Proposal , 7/76 

TURS IWT OUIPUT 
BACKOFF 




Composite, 2. dB 
Dedicated, .5dB 

WU Proposal , 7/76 

POWER ALLOCATION 

-2.5 dB 

TPM#8 

3/78 

-2.4 dB 

WU Proposal , 7/76 

CHANNEL POWER 




compos iie.e.Udc 

" WU Proposal', 7/76 - 

TORS HARDWARE LOSS 




1.2 dB 

LinCom estimate 

TORS ANTENNA GAIN 


' 


41.3 dB 

WU Proposal , 7/76 

POINTING LOSS 


-- 


0.7 dB 

WU Proposal , 7/76 

SPACE LOSS 




207.7 dB 

TDRSS Users Guide #3 

POLARIZATION LOSS 




0.1 dB 

TDRSS Users Guide #3 

ATMOSPHERIC LOSS 




0.8 dB 

TORSS Users Guide #3 

RX ANTENNA G/T 




40.3 dB |tdRSS Users Guide #3 | 
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Figure 2. 24. BER Plot For nominal Shuttle Tr .smitter. 
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igure 2.26. BER Plot for Carrier Modulator Gain Imbalance. 
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Figure 2.36. BER Plot for Subcarrier Power-Split Error. 
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Figure 2.44. BER Plot for Relative I/O Carrier Phase. 
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Figure 2.50. Sensitivity Curves. 
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Figure* 2.54.Ser.s itivity Curves for Subcarrier Parameters. 
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Figure 2.58. BER Plot for Nominal Shuttle Transmitter 
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Figure 2.63. BER Plot for Data Asymmetry. 
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Figure 2.63. BER Plot for Gain Slope. 
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Figure 2.72. BER Plot for Carrier Phase Noise (1 kHz to 150 MHz) . 
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3. LinCsim MODELING OF CARRIER AND TIMING RECOVERY 

3.1 Introduction 

This chapter presents models of carrier and timing recovery for 
the Shuttle/TDRSS S-band and Ku-band mode 1 return links In the absence 
of radio-frequency Interference (RFI ) . The generalized nonlinear link 
is described In Section 3.2. Sections 3.3 and 3.4 present expressions 
for the S-curve and r.m.s. recovered phase error for the S-band and 
Ku-band mode 1, respectively. Two sections are taken for carrier 
recovery because the two links use different loops. The carrier is 
assumed to have constant phase and frequency. Section 3.5 presents 
expressions for the S-curve and normalized r.m.s. recovered timing 
error for both NRZ and biphase (i.e., Manchester) symbols. 

The chapter following this one gives models for the RFI effects 
on the loops and gives numerical results in the form of curves for 
both no-RFI and RFI situations, assuming some typical link parameter 
values. 

3.2 Generalized Link Description 

In order to study the phase and symbol -timing recovery, we must 
characterize the signal at the input to the White Sands ground station. 
The general link diagram used herein for both S-band and Ku-band mode 1 
return links is shown in Figure 3.1. For S-band, only one symbol stream 
is present. The three symbol streams on the Ku-band return link are 
handled by treating channel 3 and the modulated square-wave subcarrier 
as one UQPSK signal and then channels 2 and 1 as another. 

Let us start by assuming that the Shuttle transmits the signal 

y^(t) * y^dj (t)cos ut - d 2 (t) sin wt (1) 
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where 

dj « unit-power stream of NRZ or biphase symbols, i * 1,2 
(It may be the output of an encoder. Nonetheless, any 
two symbols are assumed independent and are equally likely 
to have either polarity.) 

Pj = power of i— symbol stream, i = 1,2 

o) * radian carrier frequency 

Also, P.j > 0 and P 2 > 0. 

The TDR satellite filters the modulated carrier plus uplink thermal 
noise. The bandwidth of this filter is wide enough to pass the data 
undistorted; thus, the signal at the output can be written 


y 2 W s Yi(t) + n u (t) (2) 

where 

n u (t) = j/Z n ul (t)cos wt - & n u2 (t)sin wt (3) 

and n^j, n u2 are independent, identically distributed, baseband Gaussian 
processes. y 2 can be rewritten 

y 2 (t) * ^Re{[(JT 1 d ] +n ul ) + j(^d 2 +n u ^e jwt } 

= /?Re[r{t)e j(v(t)+ut) ] (4) 

where 


r 2 (t) = (» / P' 1 d 1 + n ul ) 2 + (^F 2 d 2 +n u2 ) 2 (5) 

/PZd„+n li9 vFTd,+n . 

sin v(t) = - 2 2 - - 2 , cos v(t) = ■ ( 6 ) 
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Still In the TOR satellite, the signal goes through a bandpass 
nonlinearity characterized by AM-AM and AM-PM distortion functions 
f(r) and g(r), respectively. The nonlinearity consists of an ALC/MIC 
clipper followed by a TWT amplifier. The signal at the output can be 
written 

y 3 (t) - ^'ReCf(r(t))e JCv(t)+9tr(t))+ “ t) 3 

= v7Re{[h,(r(t)) + jh 2 (r(t))]r(t)e j(v(t)+ “ t+{) } (7) 

where 

h| (r) + jhg(r) • M e i(s(lr| - 4) (8) 

The latter equation is a characterization of the nonlinearity, where 
6 is an average phaie shift introduced by g. It should be noted that 
if the nonlinearity is linear, then h 2 = 0. 

The signal at input to the ground station is 


y 4 (t) = y 3 (t) + n d (t) (9) 

where n^, the downlink thermal noise, is expressed 

n d (t) = A n d j (t)cos(wt+6 ) - A n d2 (t)sin(u)t+6 ) (10) 

The input signal has quadrature components and x 2 ; i.e., 

y 4 (t) • ^’Re{[x ) (t)+jx 2 (t)]e :i( “ t * !) } (11) 

The appropriate characterization of x^ and x 2 depends on the loop 
under consideration. 
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3.3 Carrier Tracking with Costas Loop 

3.3.1 Introduction 

The first loop we'll study is the Costas loop for use on the S- 
band return lint with its BPSK signal. 

3.3.2 Further Characterization of Link 

Since Ku-band node 1 return link will not use the Costas loop, we 
can specialize the sional to the S-band characteristics. The trans- 
mitted signal is given by (1), where now P| * 0 and d^ is a stream of 
biphase symbols. The symbol strew has two possible rates Rg, 288 Ksps 
(mode 1} and 576 Ksps (mode 2). These rates are three times the data- 
bit rates since the data bits are convolutional ly encoded at rate 1/3. 

The best way to write the quadrature components and x 2 the 
ground station input signal is as follows: 

x 1 = -(Eh 2 (r))/P^d 2 + Nj 

x 2 = (Eh 1 (r))v^d 2 + N £ (12) 

where 

*1 * -(h 2 (r)-Eh 2 (r))/PJd 2 + v ] + n dl 

N 2 * (h 1 (r)-Eh 1 (r))/^d 2 + v 2 + n d2 

v l+jv 2 = (h 1 (r)+jh 2 (r))(n ul +jn u2 ) (13) 

He see that x 2 contains most of the data and x^ just a little. Since 
most of the power entering the TORS bandpass nonlinearity is due to 
noise and not data and the noise is a much faster process than d 2 , 
then the r process defined in (5) is practically independent of d 2 . 
Similarly, h^(r)d 2 is practically independent of d 2 , i=l ,2. Thus, 
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we will assume that the noise part of is Independent of d^, 
1=1,2. However, and are not Independent. 

Ue will obtain sane statistics of the processes that will be 
useful later, namely, the correlation functions. They are given as 
follows: 




1-31^ |t |, |t| < 1/2R ? 

R 2 jt|-1, 1^ < jtj < 1/R2 

o, |t I > 1/R2 


04) 


ytj * 

K (t)P 2 + R (t) i=1,2 

n 3-i c v i c 

(15) 

11 

4J 

•f— 

06* 

(t) + R hn <*>• 1=1 » 2 

Vui Tu,3-i 

(16) 

^ ,n 2 ^ 


(17) 

R v v (t) 
V| ,v 2 

S^rVu/ 1 ^ Rh l n u2* h 2 n u2 (t ^ 

(13) 


Now that we know all about the signal that enters the loop, we 
are prepared to see what the loop does. 

3.3.3 Costas Loop Operation 

The Costas loop is shown in Figure 3.2. Input is the signal y^ 
with quadrature components and x 2 - After multiplication by 
sin(wt+0(t)) ana cos(wt+fi(t)) , where 6(t) is the loop's phase estimate, 
the signal in each arm goes through the arm filter G, which has bandwidth 
on the order of the symbol rate. When RFI is absent the power 
spectral densities (psd's) of Nj and N 2 are much wider than G, so that 
after filtering by G, and N 2 are practically Gaussian and thus not 
too hard to handle. The method of analysis of the loop is to study a 
continuous-wave (CW) phase-locked loop with the same dynamic phase error 
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process z entering the loop filter F as we have here. 

The phase error «?is defined by 

cf(t) ■ s - e(t) (19) 


It is the difference between the actual (constant) phase 6 and the 
loop estimate. Since *f varies slowly, the z process can be broken up 
into two parts: 

z(t) - S(*) + i^(t) (20) 

inhere 

$(<f) £ E[z(t)|tf] (S-curve) 

i^(t) = z(t) - S(tf), <P fixed (equivalent noise) (21) 


The equivalent loop is shown in Figure 3.3. 

S(«f) is given by 

S(tf) * \ sin(2«f)E[x|(t)-x^(t)]-cos(2«f)E[x 1 (t)x 2 (t)] (22) 

where we define for any signal a, 

a(t) = G(p)a(t) (23) 

The second term in the S-curve definition is 0 if«? has been properly 
defined so that S(0) = 0. This requirement allows us to calculate 6, 
as we will do later. The S-curve is needed for nonlinear analysis of 
the loop, which will not be pursued here. 

Our main goal is to obtain the variance of phase error. To 
do so, we assume that |«fj is almost always small so that the equivalent 
loop can be linearized. Since [1] gives full details of how this is 
done, we need only say that 
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(24) 



where H is a filter and n = n Q . 
than the bandwidth of H, we have 


Then when the psd of n is much wider 


4 ' B l NJ/(S‘(0)) 2 


(25) 


where is the one-sided noise bandwidth of H and Nq is the one-sided 
psd of n. Therefore, to calculate we need S'(0) and Ng, which we 
will later express in terms of statistics of x, and x^. 

3.3.4 Statistics of x^and iu 

Me recall (cf. (12)) that each of 3Tj and Xg consists of some of 
the filtered symbol stream plus an independent noise part, and 
Fg, respectively, the important statistics of which are given below. 


% (t) = G*(p)G(p)R. (t) (26) 

c a Z 


% * S n N (°)G*(p)G(p)«(t), i ,k = 1 ,2 (27) 

"i "i *”k 

$N i ,N k^ = J R ^^k^ t ^ dt * 1,k = 1,2 ^ 


The expression fcr (t) Is valid only when thermal noise is the only 
interference in the link. 

3.3.5 Expression ford 

Before we go on to express S * (0) and Nq, let us first find the 
angle 6 . It has been defined as that angle which allows S (0) * 0, 
which ensures that the loop has been linearized near the lock-up point. 
From (22) we find 
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S(0) * -E[xj(t)x 2 (t)] i 

so from (12), (26), (27), (28), (17), and (18), the requirement is 


seen to be that 


( 5 h,)(Eh 2 )P 2 Rj 2 (0) = 2B.S N) „ 2 (0) 


where is the one-sided noise bandwidth of the * filter. 

u 

Tnis equation is solved below for 6 in terms of statistics obtained 
by assuming 6=0. 

The first step is to note from (8) that 

hj(r) * cos(g(r)-fi) = cos 6*hj Q (r) + sin «-h 20 (r) 

h 2 (r) = sin(g(r)-6) * cos 6*h 2fl (r) - sin 6*l» l0 (r) (31) 


where 


Then 


h 1Q (r) + jh 2Q (r) = f(r)/r-exp(jg(r)> 


(Eh^)(Eh 2 ) = cos(Z5)(Eh 10 )(Eh 20 )-sin(26)^(Eh l0 ) 2 -(Eh 20 ) 2 ) (33) 

S. . (0) = cos(2s )S h h (0)-sin(26)i(S (0)-S (0)) (34) 

W h 10 ,h 20 1 "10 n 20 

S (0) = cos(2s)[S. n h n (0)-S. _ h n (0)] 

v l» v 2 h lO n ul ,h 20 n u1 h 10 u2’ h 20 u2 

- sin(26)l[S h n (0)-S. n (0)-S h (0)+S h (0)3 

^ h 10"ul h 20°ul "lO n u2 20 u2 ( 35 ) 


Writing S„ „ (0) as 
V 1 * V 2 


S„ » (0) * cos (2s) A - sin(2s )I b 
V 1 * V 2 

where the definitions of A and B are obvious, then one finds 
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tauas) > (2B 6 CA-P 2 s hi0>h20 (0n-P 2 R a2 (0)(Eh, 0 )(a 20 )}/cB G [B-P 2 (s hi0 (0) 

- ' * P 2^' 2 t0)( ^^Q^* ( ^h 2G ) ^ ) ^ < 37 > 

3.3.6 Expressions for S'(0) and Nq 

Now that we have the 6 that properly defines cp, we may proceed to 

obtain the quantities S'(0) and needed In the calculation of the 

2 

phase error variance o«p given in (25). Since now 

S(«f) = \ sin(2V)E[x|(t)-x^(t)] (38) 


then the derivative at 0 is 

S'(0) = R- (0) - R- (0) 

x 2 x, 

s ((Eh 1 ) 2 -(Eh 2 ) 2 )P 2 R a2 (0) + (S^(0)-S M (0))2Bg (39) 

The equivalent noise n(t) and its two-sided psd Nq/ 2 are given by 
n(t) = -x 1 (t)x 2 (t) 

= (Eh 2 )^ff 2 (t)?T 2 (t) - (Ehj )/P^d 2 (t)lt| (t) 

-^(t)N 2 (t) + (Eh 1 )(Eh 2 )P 2 3 2 (t) (40) 

N 0 /2 £ p 2 R a 2 (0 ) [ t Eh l )2s N 1 (0)+(Eh 2 )2s N 2 (0, ‘ 2(Eh l ){Eh 2 )S N 1 ,N 2 (0):i 

+ % W ^ ) 

" 1 W 2 


Since in the absence of P.FI Rj and N 2 are usually not very correlated 
and since the G filter is RC-type, it is sufficient to approximate 
% w (0) by 

V2 
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Vz ' 0 ’ 4 B <y o) y o) <«> 

3.4 Carrier Tracking with Two-Channel Costas Loop 
3.4.1 Introduction 

The Shuttle Ku-band mode 1 return link will use a two-channel 
(i.e., quadriphase) Costas loop with hard-limiters to track each of 
the carrier and subcarrier. In order to study both cases with the 
same model, we keep the full generality of the link model shown in 
Figure 3.1 and of the signal characterization given in Section 3.2. 
3.4.2. Further Characterization of Link 

Since now both symbol streams dj and d^ are to be used for carrier 
tracking, the most useful expressions for the quadrature components x^ 
and Xg of the loop input signal are different from before; they are 

X| * (Ehj)i/Pjd| + N-j 

x 2 = (Eh 1 )/P^d 2 + N 2 

* (h-|(r)-Ehj)^yl.| - h 2 (r)^d 2 + v 1 + n rfl 

N 2 = (h^rJ-Eh^ )yf^d 2 ♦ h 2 (r)v / P r Jd 1 + v ? + n d2 (43) 

where v.,h., i = 1,2, are as defined before in (8) and (13) and now£ 
is defined by 

« = tan _1 [E(f(r)/r-sin g(r))/E(f(r)/r-cos g(r))] (44) 

so that 

Eh 2 (r) = 0 (45) 

It will be shown later that this definition of 6 is adequate. 
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Just as for the Costas loop, we assume that the symbol streams 
dj and dg are Independent of Nj and Ng. For channels 2 and 1 (those 
on the subcarrier) this is practically true by the same reasoning as 
for the Costas loop used at S-band. This applies also to roughly half 
of the transmission on channel 3 (on the carrier), the most powerful 
channel. For the highest data rates or for the highest values of 
E b /N 0 u (data-bit energy to one-sided psd of uplink noise) which will be 
used for any particular data rate on channel 3, we must appeal to the 
argument that hj(r)-Ehj and h 2 (r) are small in order to justify the 
assumption. 

We will need later the correlation functions of d. and N. , given 
below. 


K (t) = 



{t | < l/R. 

o. 

it! > l/R. 

l-3R i |t|. 

|t| < 1/2R. 


1/2R. 1 |t| < 

o , 

|t|>l/R. 


(NRZ) 


(46) 


N 

R N (t) * R h .(t)P 2 + R h _(t)P 1 + R v (t) +-0d 6 (t), 1 = 1,2 

(47) 


R N r N 2 (t) * (P 1 


P 2 )R h,.h 2 


<t) * R, , (t) 

V V 2 


(48) 


where R^ is the symbol rate of d^ , i = 1,2, and where correlation 
functions of v^ and v 2 are the same as given in (16) and (18). 

We now know what the signal is that enters the loop, so we can 
start to look at the loop. 
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3.4.3 Operation of Two-Channel Costas Loop 

The two-channel (l.e., quadriphase) Costas loop with hard-limiters 
that will be in the ground station is shown in Figure 3.4. As opposed 
to a more optimal implementation [1], this implementation uses the 
same filter 6 in all the arms. The ratio Cg/Cj will be 2 when P 2 /Pj = 4. 
Just as for the Costas loop, the method of analysis is to study an 
equivalent CW phase-locked loop with the same dynamic phase error 
process e input to the loop filter F. 

Since the operation of this loop is fully described in [1], we 
will go over it only briefly. The phase error*? is defined as in (19). 
The e process is the sum of two others, as shown in Figure 3.4: 

e = -ej + e 2 (49) 

Hie S-curve S(«?) and the equivalent noise n^(t) are given by the 
following: 

S(«f) = E(e|«?) (50) 

n^(t) = e (t) - S(ef) for fixed 6 ? (51) 

We obtain below the S-curve and the quantities S * (0) and Nq for 
the equivalent loop. The latter two are used to calculate the variance 
of fusing (25). All are expressed in terms of statistics of Xj and x 2 . 

3.4.4 Properties of x^ and x^ 

The bandwidth of the G filter is chosen to pass most of the 
power of both d^ and dg. Unlike for the Costas loop, we need to assume 
here that the effect of G on and on d 2 is merely to reduce their 
power; i.e.. 
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x^t) * Ajdj + TTj 
A i * ^(Eh-jJi/F^ 

n 1 * Rj^(O) (52) 

In the absence of RFI R^ and fig are approximately Gaussian. The 

correlation functions of ffj and R^ are as given in Section 3.3.4. 

3.4.5 S-Curve Sfo) 

Now we are able to write an expression for S(«(>), which from (49) 
is seen to be the sum of two terms. 

S(«?) = + S 2 («f) 

Sjfo) = E[e.j(t)|*f] (53) 

From [1] we obtain 

E l^l * [“Sin^-tAi dj +Rj l-cosY-fA^ d 2 +TT )]sgn[cosY-(Aj dj+Rj ) 

- sin'MAg <» 2 +^ )] 

= [-sinYAj dj-cosY/^ d 2 +(-sinY*Rj -cosVf^ )] 

• sgn[cos«PAj dj-sinY-A ? d 2 +( cosYRj -sinY-^ )] ( 54 ) 

so that for correlated Rj and R 2 , 


s 1 («e)/c 1 


- ^(sinY-Aj 

- ^(sinY-Aj 


+cosY- /^)erf 
-cosY-/^ )erf 


( 

( 


cosY-A^ -sinf-Ag 


J2xs-y 

cosY-Aj +sinY-A 2 

a-, 


) 

) 
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J.s1n(*)(oj^ -o|, ) - cos(2*)Rj^ jj^(O) 

* 

-(cos^-A. -sin^-A, ) 2 /2o? -(costf-A, +sin<PA 2 
-[e 1 + e 


) 2 /2a 

(J 


where 

a* * cos 2 (cf)a^ + sin 2 (^)o^ - sin(2*)R^ ^(0) 

and the facts have been used that for Gaussiar. random variables u,v 
and constant A, 

E sgn(A+u) = 

E[v sgn(A+u)] 

erf(x) s e' y dy (E 

J 0 

Similarly, e 2 is given by 

c 2 /C 2 = [cos £ P*(A 1 d 1 +H 1 )-sin«f.(A 2 )]sgn[-sin < P.(A 1 d^ ) 

-cos‘P-(A 2 d ? +F 2 )] (E 



so that 

SjjW/C., 


/sin^-A, +cos‘?A 2 

- ^cos'P-A 1 -sin^-Ap )erf ^ ^ 

/sinV A, -cos^-Ap 

- ^cosf-ft, +sin<f-A z )crff- jr^r 


) 

) 
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C-j A- j+C^A v and S(0) is much less. 

3.4.6 Psd of Equivalent Noise 

Still needed to evaluate (25) for the phase error variance is 
the one-sided psd Nq of n(t) s n Q (t). One finds 

n(t) * -e-j + r.g + S^O) - S 2 (0) 


= ^(A^+N^sgnfA^+flj) - )sgn(A 2 d 2 +ff 2 ) 

+ S 1 (0)-S 2 (0) (64) 

OD CX> 00 

NA/2 = f R (t)dt + f R (t)dt - 2 f R (t)dt (65) 
u j e, J e 2 ' e l* e 2 

•00 .00 .OO 

The expression for Nq in the case of ^ ,ffi 2 uncorrelated is given in 
(79) of [1]. 

3.5 Symbol Tracking with Digital Data-Transition Tracking Loop 
3.5.1 introduction 

The digital data-transition tracking loop (DTTL) will be used 
for both NRZ and biphase symbol synchronization on the Shuttle/TDRSS 
S-band and Ku-band mode 1 return li , The loop input is one of the 
quadrature components and x 2 of . .e ground station input y^, given 
by (12) for S-band and by (43) for Ku-band mode 1. Figure 3.5 shows 
the loop, the exact operation of which is described in [2]. The input 
to the loop is of the form Ad(t+e)*w(t) , where d is a unit-power 
stream of NRZ or biphase symbols and w is white Gaussian noise of 
one-sided power spectral density (psd) Nq. We obtain below the S-curve 
and the noise psd tor an equivalent CW phase-locked loop, shown in 
Figure 3.3; then we compute the sigma of the normalized timing error 
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for the linearized equivalent loop. 

3.5.2 Tracking NRZ Symbols 

The performance of the DTTL on NRZ symbols is extensively analyzed 
in [3j. The S-curve and equivalent-noise psd expressions given below 
come from that source. 

The notmalized timing error A is defined by 


A = (c-2)/T (56) 

where A is the timing of the data signal, e is the timing estimate 
produced in the loop, and T is the symbol time. The S-curve g( x) 
of the equivalent loop is defined by 

g(A) = E[e k |A] (67) 

and is given by 



. A erf(^(l-2A)) - ^ U 0 -2 a] 

• [erf(4i^) - erf(^(l-2A))] for A *4p 

^ erf(^(l-2A)) for-^<A<^ 

( 68 ) 
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The two-sided psd of the equivalent noise when x = 0 is given 
J>y 



- 7? ^ «p(-R s )erf(^) - ^ R s erf 2 ^) {71 , 
Frcmi (68) , we see that 

9^(0) = erf(/BJ) - ^ exp(-R s ) (72) 

t 

so we find that for the linearized loop 



where B L is the one-sided loop bandwidth. In order to compare 
later the loop's performance for NRZ and biphase symbols, we note 
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3 -5.3 Tracking Bi phase Symbols 

Me assume that tfie UTTL tracks a stream of biphase symbols as 
If It were a stream of double-rate NRZ symbols. Figure 3.5 still shows 
the loop operation If T is replaced by T/2 and If now Cg Is a fraction 
of T/2, the duration of half a biphase symbol. 

The normalized timing error X is still defined as in (66). The 
S-curve g(x) now has period 1/2. Starting from Equation (A-14) in pj 
we find the S-curve: 



Its derivative at zero is given by 


9 o<°) 




(76) 


From Equations (B-37) and (B-42) in p ] after some work we obtain 


E[ e k e k+ffl !x=°] * 0 for all m > 1 (77) 

so that the psd of the equivalent noise for X * 0 is given by 

S(0,0) = j E[e k Jx=0) (78) 

From Equation (B-30) in [3] we obtain 
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where R Is defined as in (69). Then o is given by (73). 

S A 

In the case of low noise, for the same values of £q» R $ , and B^T 
biphase symbols will give a smaller timing sigma than NRZ symbols will, 
as can be seen from (74) and the fact that for biphase symbols 


lim 

v~ 


w i 

tfK' 7 * 


(80) 
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4. RFI EFFECTS ON SHUTTLE/TDRSS LINK 

4.1 Introduction 

The Russian air defense system uses large numbers of powerful 
radars whose carrier frequencies coincide with the TDRSS SSA return 
link frequency band. The TDRSS payload and ground station hardware 
had to be redesigned to minimize the deleterious effect of the radio 
frequency interference (RFI) resulting from these radars. The overall 
effect of the RFI and the hardware changes is impaired system perform* 
ance both with RFI and without. 

This chapter gives a statistical description of the RFI environment 
and lists the hardware changes in the TDR satellites and ground station. 
Then analytical models are described to assess the performance impact 
both for the bit error probability and for carrier and data tracking. 
Finally, preliminary performance results are given. 

4.2 RFI Environment 

The radio frequency interference is caused by large numbers of 
high-powered radars (peak EIRP in excess of 100 dBW) in eastern Europe 
and Russia. The radar region as seen from the orbital locations of 
the two active TDRSS satellites is shown in Figs. 4.1 and 4.2. The 
illumination of the TORS occurs through the radars' main, side, and 
backlobes, resulting in a large dynamic range of the composite RFI. 

The radar pulses consist of pulsed sine waves (CW) accompanied 
by wideband Gaussian noise (Fig. 4.3). The pulse duration is in the 
2 to 5 microsecond range. Table 4.1 shows an unclassified coarse 
approximation to the actual distribution of pulse power and pulse 
duration. The first two lines correspond to radar pulses whose CW 
component lies outside the TDRS input band; hence, only the wideband 
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noise affects performance. The remaining lines represent CM pulses 
within the TORS band. Since the CU power by far exceeds the wideband 
noise power in this case, the latter my be neglected. The first 
column shows the RFI EIRP in the direction of the TORS (for noise-like 
pulses it is measured in a 20 MHz bandwidth) quantized in ten-dB steps. 
The pulse durat on has been quantized into two values, two and five 
microseconds. The numbers in the second and third column represent 
the number of pulses per second with a given power and duration. The 
duty cycle listed in the fourth column is the product of repetition 
rate (PPS) and pulse duration. The environment shown in Table 4.1 
applies to the TORS over the Pacific and for the center frequency 
2217.5 MHz. It is generally referred to as the non-benign Shuttle 
environment [1]. The RFI environment encountered by the eastern 
TORS is worse, the other Shuttle S-band frequency, 2287.5 Wz, is 
less severely affected on both satellites. Table 4.2 shows a more 
severe TDRSS SSA RFI environment which is being used for program 
testing. 

Hie TORS S-band antenna pattern can reduce the interference 
when the Shuttle is not too close to the RFI region. The antenna gain 
drops about 12 dB at 1.5 degrees off-pointing and about 24 dB at 
4 degrees off -pointing (Fig. 4.4). Figs. 4.1 and 4.2 show where 
those off -pointing angles are reached. 

4.3 TDRSS RFI Hardware Changes 

The TDRSS RFI hardware modifications have the purpose of removing 
powerful spikes from the signal. This is accomplished by clipping the 
signal when it exceeds a certain level. 
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Figure 4.2. TDRS East View of Radar Interference Region. 
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Table 4.1. Moderate SSA RFI Environment. 


Interference 


S IMPLIED TEST ENVIRONMENT 
(PACIFIC TDRS) 


Moderate SSA (Including Non-Beni gn Shuttle) 
PPS Duty Cycle 


dBH 

5 jusec 

2 jusec 

(percent) 


25 

2000 

8000 

2.6 

Noise 

35 

6000 

2000 

3.4 

Noise 

45 

2000 

300 

1.1 

CW 

55 

2300 

100 

1.2 

CW 

65 

250 

50 

.14 

CW 

75 

25" 

5 

.001 

CW 

no 

1 


N/A 

CW 


- 130 - 


■JbnC* 




jCinCi 


om 


In the satellite this clipping is performed by a limiting micro- 
wave Integrated circuit (MIC) amplifier located in the automatic level 
control (ALC) circuit. This set-up is shown in Fig. 4.5. The 
characteristic of the MIC amplifier (Fig. 4.6) shows good linearity 
below the clipping level and a sharp cutoff at this point. Presently 
the nominal operating level is set 6 dB below saturation but this 
setting can be varied between 2 and 10 dB below saturation in steps 
of .5 dB. 

Unfortunately, the limiter also introduces nonlinear phase 
distortion in the form of AM/PM conversion. This characteristic is 
shown in Fig. 4.7. 

There will also be a clipper in the TDRSS ground station. However, 
its characteristics and nominal operating level are not defined at 
present. 

4.4 LinCsim Modeling of RFI Effects on Bit Error Rate 
4.4.1 Introduction 

This section presents the models and algorithms used in LinCsim 
to compute the bit error probability of the nonlinear TDRS channel 
in the presence of RFI. 

The RFI environment as presently characterized for the TDRSS 
consists of RF pulses, either noise-like or pulsed CW, of approximately 
2 to 5 psec duration a.id a wide range of power levels. Typically, a 
large number of the pulses exceed the signal power level. The duty 
cycle of the RFI, i.e., the product of pulse duration and number of 
pulses per second, is used to characterize the severity of the 
disturbance. Typical values are between 20 and 30%. 
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Figure 4.7. ALC/MJC AM to PM Conversion Characteristic. 



The basic approach is to condition first the error probability 
on the satellite repeater output. (The conditioning is not done on the 
complete repeater output waveform but rather on a parameter which 
provides a sufficient statistic for the signal in the detection process. 
This parameter, called the decision variable Zj M below, is the demodulated 
repeater output integrated over a symbol time.) Since the only random 
disturbance left is the downlink noise this conditional error rate is 
given by a simple expression involving the standard error function. 

This conditional error rate is then averaged over the statistics of 
the decision variable. 

For this averaging of the conditional error rate different 
approaches are used for low and high data rates, where the break- 
point is the inverse of the RFI pulse duration. The first approach 
is to compute the probability density function (p.d.f.) of the 
decision variable using the classical moment technique and then 
to average the conditional BER over this p.d.f. This applies to 
the case of low data rates since the RFI pulse duration is only a 
fraction of the data symbol time. The second approach, used for 
the case of high data rates, is based on the assumption that one 
bit is only affected by at most one RFI pulse. The error rate is 
therefore computed conditioned on the presence of a particular RFI 
pulse characteristic and then averaged over the probability distribution 
of the pulse characteristics. 

The detailed description of these two approaches is given in 
Section 3. 4. 3,4. 4. 4. Section 4 A2 defines the models for the channel 
and for the RFI and defines the notation used in the remainder of the 
chapter. Section 4A5 addresses the problem of computing the moments 
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needed In the approximation of the decision variable p.d.f. These 
moments are derived from the characteristic function, in part 
analytically, in part through fast ntmierical algorithms, depending 
on the character of the random variable. 


4.4.2 Description of Model 
4.4.2. 1 RFI Model 

The RFI environment as seen by the TDRSS satellites consists of 
two classes: pulsed wideband Gaussian noise and pulsed CW tones. 

These two classes are further divided into groups with different 
power levels. The pulses from each such group are assigned to arrive 
at the TORS as a Poisson process and independently from all other 
groups. He will however assume that no pulse overlaps occur which 
is true with good accuracy for duty cycles up to 30Z and provides 
an upper bound on the error rate for higher duty cycles since 
effectively a larger portion of each symbol is affected by RFI. 

Hie pulsed RFI is represented as follows 


where 


j(tj - J c (t) ♦ J 9 (t) 

- "c n K 

: (t) * E E J k r )P(t-t k , ) 

k =1 *1 c c c c c 


yt> • 


c ' ~c 
»- "k. 


E E J C (t,p(t " t k v )exp(J “o t) 
5 90 


9 9 

n(t) * a complex baseband Gaussian random process 
with mean zero and variance 1/2. 


0 ) 

(2a) 

(2b) 


r k (t) 


1 0 < t < i. 

0 elsewhere 


(3) 
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x 8 Pulse duration 


„ /number of groups of the RFI for CM tone and noise-like, 

"g / trespectively. 

t 

n kc) /number of RFI pulses in group k (or k 1 arrived in 
ft t /= \the observation interval (tg,t)9 9 

k 9 

t 0 = the beginning of the observation time, usually the 
beginning of a symbol. 

tj i * the arrival time of the ith RFI pulse in the ith group. 

Here, the variables n^ and n^ are homogeneous Poisson random variables. 

The CM tone RFI will be further expressed as follows: 

H c . n kc 

J c<‘> 1 E E J k (tt k t ,exp( *k t 1 exp ( J V) 

k =i t =1 = “ * * 


Jet*? * 


c c 


= jJJ exp j(wgt+*j) 


where 

2 /"c ""c 

i J J ■ E E 

V' V 1 


* > tos *k, 

c c c 


/ N c "kc V 

E E J k p(t - t k . > sin ♦* « 

V => c c c c 7 

c c * 


\ "kc 


.(t) _ ...-1 *c ' *c 
*J ‘ ta " N c n k( . 


E E V (t -v » >*<« 'u r 

r*i r^i c c c c c 


£ £ Vc 


c c 


Mote that |J C | and 4>,(t) are functions of time because of the pulse 

p(t-t k t ). However the variations of j J | and ^j(t) are much slower 
c c 

than the carrier frequency and they remain constant over the 

RFI pulse duration t. 

cllinCoi 
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4.4. 2. 2 Channel Model 

The mathematical model of the TORS link is depicted in Fig. 4.8 

where there is RFI in the uplink channel and Gaussian noise in the 

uplink and downlink channels. The signal transmitted is a BPSK 

signal or a QPSK signal modulated with independent inphase and quadrature 
data streams, so at the receiver the data messages are detected 

for the inphase and quadrature phase separately. In this study 

we shall discuss the BPSK signal case. The signal received at the 

satellite transponder can be expressed as 


x(t) = V exp(j(o. 0 t+e(t)) + J(t) + n u (t) (5) 

where 

o 

v / 2 - signal power 
e(t) = transmitted phase 
ojq/2u = carrier frequency 
J(t) = RFI, as defined by eq. (1) 

' n u (t) = uplink channel noise with mean zero and variance a* 

Using the RFI model of eq. (1), the expression for x(t) is rewritten 
as 


x(t) = V exp(j(« 0 t+9(t))) + J c (t) + nj(t) (6) 


where 

n l(t) * n u (t) + J g (t) 


Note that the noise process n^(t) still has zero mean but its variance 
changes with time (conditioned on the RFI arrival process) because 
of the presence of the noise-like RFI. By the assumption that the 
channel noise process and the RFI are statistically independent of 
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each other, the variance of n^ft) may be either equal to o or 

2 2 U 
® u * Jjj /2, depending on which group of RFI is present. Here we assume 

that samples taken from these noise processes at the Ryquist rate 

are statistically independent of each other with variance o 2 if no 

2 2 u 
RFI is present, or o u + **kg/2 ^ RPI ^gt^ group present. 

Now the signal x(t) passes through the satellite transponder 

and becomes 


y(t) = f(R)exp(j(a. 0 t+e(t) + n + g(R))) (7) 


where 


R * the envelope of x(t) 

n = the phase of x(t) 
f(R) = AM/AM distortion 

g(R) = AM/PM distortion 

One can show, [see 2] that the probability density function of R and 
n, after averaging over the random phase of RFI, 


P(R,n,|J c |) = P 1 (R,n)f(R,n,(J c I) 


where 


P^R.n) = 


2vo 

f(R»n»|J r |) = exp 


exp - -U (R 2 +V 2 -2RV cos r) 

2moj 2<Xj 


( 8 ) 


(9a) 


(■ M 

1 f* 

Iq(x) ■ - I exp(x cos e)de 
Jn 


/r 2 +V 2 -2RV cos 


1 

(9b) 


If |J C | equals 0, the expression P(R,n,0) reduces to that for normal 
channel condition, that is 


P(R,n,0) * P-j(R,n) 


00) 


~<jCinC^om 
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4*4 .3 Conditional Error Probability 

At the ground station, the received signal can be represented 

as 

r(t) * Re[y(t)] + n<| (t) 

where n d (t) is the downlink Gaussian noise with mean zero and 

o 

variance in the receiver channel band. The signal is first 
demodulated by a local carrier 2 cos(u» 0 t+g), ideally tracking the 
received phase of the signal and processed through an integrate-and- 
dump (I&D) device which gives 

Z. = J- f d f(R)cos(e+n+glR)-g)dt + n (11) 

X T d J o • 

where 

n z '= baseband noise process whose mean is zero and variance 

4 is <°d/ T d> 

= integration interval of the I&D, assumed to be the same as 
the data symbol interval. 


The integral term represents the desired signal for the detection in 
the absence of uplink interference such as channel Gaussian noise 
and RFI. Denote the power associated with this term as E d , then the 
downlink bit energy to r.oise spectral density ratio (E^/Nq) is defined 
as where o £ is the variance of the downlink noise in the data 
bandwidth. 

When the uplink interference is present, the power E d is shared 
by the retransmitted signal and interference. Therefore, the effect 
of the uplink interference is two-fold— power robbing and signal 
perturbation. Since the integrand in (11) is a nonlinear function 
of the uplink signal, noise and RFI the statistics of the associated 
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integral cannot be found exactly. However, this integral can be 
well approximated by the sum of signal samples taken at the Nyquist 
rate of the repeater input signal : 


where 



+ n_ 


02 ) 



(13) 


N = product of the channel bandwidth Bjp prior to the 
nonlinearity and the data symbol time 


a t = f(R { )cos(n 1 -»g(R i )-i) 


(14) 


The samples a., taken every B“p sec, are statistically independent 
of each other, but not necessarily identically distributed because 
of the possible presence of RFI. The impact of an RFI pulsed on 
the variable Zj^ can be classified into two cases— (1) all a^ are 
affected. by an RFI pulse, (2) parts of the samples of Zj M are affected. 
Oenote the pulse duration of an RFI pulse by T . If Bjpt is smaller 
than N, i.e., the data symbol is longer than the pulse duration T , 
only some of the samples a^ are affected by the RFI. We shall group 
those independent and identically distributed ( i id ) a^ into groups 
as follows 


- .1 

IN N 


N ft, 


E E 

k -1 rSl 


a r tfc c> + 


\ % n 0 

£ E MV ♦ EM 01 

r=l / 


V 1 r ' 


5) 


where 
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For a high data rate, all of the samples are completely hit by a 
single RFI pulse. So the variable Zj M for symbol detection is 

z .h ■ jr £ MV ‘ 16 > 

1=1 

conditioned on the k th group of RFI pulses. Equations (15) and (16) 
point out the difference between the low and high data rate models 
alluded to in the introduction. In both cases the error probability 
conditioned on the decision variable Zj^ takes the form 

W - -5erfc(^-) (17) 

•where 

erfc(x) f exp(-t 2 )dt 
7 •'x 

The averaging over Zj M however is done differently in these two cases. 
This will be explored in the next section. 

4.4.4 Bit Error Probability 

Due to the different impact of RFI pulses on the decision variable 
Z JM , two approaches are needed to evaluate the bit error probability 
for high and low data rates, respectively. 
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4. 4. 4.1 Low Data Rate Case 

In principle, the error probability for low data rates can be 
obtained by averaging the conditional error probability oi' eq. (17) 
over the statistics of the variable Z JM , which is not explicitly 


known. 


p e * E Z 


IM 


0.5 erfc 




( 18 ) 


However, we shall use the classical moment technique to construct 
the approximate probability density of Zj M - To do this we need to 
evaluate the moments of Zj M with respect to the statistics of the RFI 
pulses and their arrival distribution and the uplink channel noise. 
For example, the kth moment of Zj M can be found in a straightforward 


manner as below: 
m. = E 


'./"c V "c\ n 0 

7 E E* V* EZ?V + E 

" \k c =l r=l c k =1 r=l 9 r-1 



(19) 


One can recognize that the evaluation of the above equation is non- 
trivial. Thus the characteristic function will be used to evaluate 
the needed moments, as discussed in Section 4.4.5. Based on the 
calculated moments of the classical moment technique, the expectation 
of the error probability in eq. (18) can be written as 

v / Z 1 

'2 


P E * 0,5 2 


erfc 


k=l 


(tsst) 


( 20 ) 


where I s a discrete probability density function (p.d.f.) 

approximating the continuous p.d.f. of the variable Zj f/) and v is 
the number of points Z lk with nonzero probability 
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4.4.4. 2 High Data Rate Case 

Since the duration of an RFI pulse becomes longer than the 
symbol duration for a high data rate, the error probability 
conditioned on a type of RFI pulse should be computed, then averaged 
over the arrival distribution of the RFI pulses. Denotp by P (k ) 

C Xr 

the bit error probability conditioned on k^-th group of RFI pulses 
Then the overall bit error probability is (see Appendix) 

N c N g 

p e = £ P(k c )P e (k c ) ♦ £ P(k g )p e (k g ) ♦ p 0 p 6 (o) (2i) 

V 


where 


v 1 


:!!■ 


p(k c ) L 

?(k. 


Probability of a symbol being hit by 
pulse of k c ~th (or k^-th) group 


an RFI 


Pq = Probability of a symbol not hit by any RFI pulse 

P (0) = Error probability conditioned on the absence of RFI 
e pulses 


One can express Pq as below: 


- t p(k c> + E p(k 9>i 

W 1 k - =1 


c 9 

The probability P(k ) takes the form 


(22) 


P(k t ) = A k T k (23) 

x i 


where 


= Arrival rate of the k £ -th group RFI pulses 
i 

t. = Pulse duration of the k th, group RFI pulses 
Here we have assumed that the overlaps between RFI pulses are negligible 


and that X. is small (less than one), 


cjCinC^otn 
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4.4.5 Evaluation of Moments 

In computing the conditional error probability for high data rates 
or overall error probability for low data rates, the moments of the 
decision variable Zj H are needed in constructing the probability density 
function of Z JM . In this study, the characteristic function approach 
is used. 

Denote the characteristic function to be * 7 (jw) 

*IM 

♦ z (j«) = E[exp(j w Z IH )] (24) 

IM 

■ £ <»» 

k=0 

where 

E = expectation taken over all randan variables 
ra^ = k— moment of Z JM 


A comparison between eq. (15) and (16) shows that the high-rate 
model can be evaluated as a simple special case of the low-rate 
model, hence only the averaging for the low-rate case is discussed 
below. 

Substituting eq. (15) into eq.(24) and based on the assumption 
that samples in each group are statistically independent and 
samples among one group are iid, we can write the characteristic 
function as follows: 

N. N 

♦ z (j») • *!J(j g) TT \ (i p "IT \ (j t> ( 26 ) 

L W u N k c *l k c K k g =l K g ” 

where 

* 0 U*>) * E[exp(j w a r (0)] (27) 
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sV J "> * ♦srfcr ECexp0 “ a r (k 9 ,,] 

\ Ua) ' E js\ < J “>j 

I \ I 

♦ k (j«) s Ej s * k 9 (je)| 


(28a) 


(28b) 


(29a) 


(29b) 


Here we denote ^ (j u ) (or (j«)) to be the characteristic function 
c 3 

of a single sample a r (k c ) (or a r (k g )), divided by the characteristic 

function of a single sample a r (0), and (ju) (or ( u )) to be the 

c g 

resultant characteristic function averaged over the statistics of 
the Poisson arrival rate of the* k c -th. or (k^-th) group of RFI pulses. . 
For a given data rate R, the total number N‘of samples a (k ), 

I Ju 

k = 0,1,2,..., can be written as 

Mr 

N = B if /R (30) 

Thus the variable n. , defined by Eqs. (15a) and (15b) can be 

K £ 

expressed as 


"k. - 


and Eqs. (29) become: 


( f NRx 

\ ( J“> ■ E |[s\ P 'J“> 


1 ° k I 

] c l 


*k (j «) * E j[s*k 


(32a) 


(32b) 


Here n£ and n£ are the only random variables and their probability 
c g 

distributions are characterized as Poisson processes: 
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K (n) 


(PJ)" “ P . T 

nl e 


where 


P t = Pulse repetition rate of ith group of RFI 
T * Observation period. 

To simplify the notations, we denote 

NRt 

H k (ju) 5 s « k P (j jf) 


H 0 (j») - «"(j |) (35) 

Then, after taking the expectation with respect to nT for Eq. (32a) 

T c 

or to N£ for Eq. (32b), we have 

9 

*k (j V " exp(A k T d (H k (*) 

1 £ £ 

In general , the functions H. (j«) can be expanded in terms of power 

£ 

series in (ju) and HqUu) in terms of its semi -invariants: 

(j«) - ]£ m^kj) (37a] 

1 i=0 


H 0 (ju») = exp 




Note that the variable m^k^) is not the moment of some physical 
random variable. However it is a function of the moments of the 
sample a (0) and the moments of the sum of the samples a (k ), 

\ \ 

either £ a r (k c ) or V a r ( k g ) 

r=l r=1 

Now using Eq. (37a) and (37b) we have 


cHux(^o\ 
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* z IM - exp £ (VO) * £ P k T /i( k c ) 
1*1 k c »l c 


♦ £ VA'V 1 


k = 1 9 

9 

This yields the resultant i-th semi -invariant of Z, 


X, = VO) * £ p.yvv ♦ £ P k Vide) (39) 


k =1 c k^1 9 

c g 

for i = l,2,3,...m 

Therefore, the moments of Z JM can be recursively computed from the set 
of semi -invariants Xj as follows 

\ = k 1 

i=l 


_ ■ n! 

V HTn^iTT 

4.5 LjnCsjm Modeling of RFI Effects on Synchronization Loops 
4.5.1 Introduction 

RFI will degrade carrier phase and symbol timing tracking when 
it is present on the Shuttle/TDRSS S-band return link. In this section 
the RFI environment of Table 4.1 is assumed. The duty cycles are so 
low that we can well assume that no pulses overlap. The way in which 
the RFI will affect each tracking loop depends on the relationships 
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among the symbol rate, RFI pulse durations, and RFI pulse repetition 
rates, tfe will give new expressions for r.m.s. phase error and timing 
error, to replace those presented in Chapter 3, that take RFI into 
account. 

4.5.2 RFI Effect on Ground-Station Input Signal 

Let us examine what the RFI effect is on the signal at input to the 
ground station. We refer to the generalized link shown in Figure 3.1. 
Since the correlation time of the TORS input filter (about .05 usee) is 
much less than the pulse durations, that filter passes the RFI essentially 
undistorted except for bandlimiting it. The nonlinearity (combination 
clipper and TUT) is memoryless. Thus, at the input to the ground 
station a set of signal statistics something like those in Section 3.3.2 
can be found conditioned on each RFI situation. Since the power levels 
of the WGH RFI are relatively low, during no RFI or a WGN-RFI pulse 
the clipper doesn't have much effect. The noise power entering the 
TUT is just somewhat larger during a WGn-RFl pulse. However, the power 
levels of the CW RFI are much higher than those of the data and uplink 
noise, so that the clipper suppresses the data; little more than a tone 
enters the TWT. 

We now proceed to treat the carrier and symbol tracking loops 
separately. 

4.5.3 RFI Effect on Carrier Tracking 

4. 5. 3.1 General Description of RFI in the Loop 

We need to model the effect on the Costas loop, shown in Figure 3.2, 
of each RFI process. Fortunately, the general behavior is the same for 
all since the pulse duration and repetition rates of all bear the same 
relationships to the two S-band symbol rates. The meaning of this will 
be made clear below. 
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First step in the modeling is to consider the response of the 
arm filters 6 to a signal with a RFI pulse. The one-sided noise 
bandwidth Bg of G is .471 Wz for S-band mode 1 (symbol rate 288 Ksps) 
and .942 Wz for mode 2 (syntool rate 576 Ksps) [3], corresponding to 
a correlation time of l/4Bg, which is '.53 psec and .27 nsec for 
modes 1 and 2, respectively. This is much less than the pulse duration, 
especially since most of the higher-power RFI pulses, the ones that 
matter the most, have duration 5 psec. Therefore, we assume that at 
the output of each of the arm filters, signal statistics can be obtained 
conditioned on each RFI situation. 

The next step is to realize that since the loop merely multiplies 
the outputs of the arm filters to obtain the dynamic (base error process 
z, then the statistics of z are similarly conditioned. 

Therefore, because of the definition (3-21) of the S-curve S(«f), 

S' (0) for use in the calculation (3-25) of phase-error variance is 
given by 


S*(0) = [S' (0) Jno RFI]-Pr[no RFI] 


+ [S* (0) | in a pulse of kth WGN-RFI process] 

k=l 

-(duty cycle of kth WGN-RFI process) 

"c 

+ [S'(0){in a pulse of kth CW-RFI process] 

i&i 

.(duty cycle of kth CW-RFI process) (40) 
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where it is assumed that the pulses of each RFI process occur at a 
constant rate and where My and are the nun.be rs of UGN- and CW-RFI 
processes, respectively. Since phase error V is just scaled 2 
passed through a filter H of very long correlation time (cf. (3-24)), 
then Nq for use in (3-25) is given by a similar equation to (40). 

He must now obtain conditioned quantities. 

4. 5.3.2 S'(0) and Nq During Wo RFI or WGN RFI 

The situation during a WGN-RFI pulse is like that of no RFI since 

WGN RFI looks like increased uplink thermal noise. Therefore, S'(0) 

and Nq conditioned on one of these situations are calculated from 

(3-39) and (3-41) using the conditioned values of Eh. and S N (0), 

i =1,2, of S u at* (0) and of S^- » (0). Since the clipper doesn't 
n*| ”2 

do much during HGfl RFI, then Nj and N 2 are still almost independent, 

so (3-42) can be used to approximate Sjj » (0). 

N 1 "2 

4.5. 3. 3 S'(0) and Nq During CW RFI 


The signal arriving at the ground station during a CW pulse is 
essentially just a tone of constant power. This is because the RFI 
power is at least 15.7 dB above that of uplink thermal noise in the 
bandwidth of the TDRS input filter and the noise power is greater than 
that of the data. The TDRS clipper cuts down the power to a fixed level 
and the TWT amplifies the signal. The downlink noise is inconsequential. 

Let 'j*: consider what the loop's arm filter 6 does to the tones. 

Since the correlation time of G is much less than the pulse durations, 
the tones appear to G to have infinite duration. The filtered quadrature 
components and Xg of the loop input signal during the ith^ pulse are 
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X|(t) = G(p)x 1 (t) = /CTTpcosfu^^.) + n dl (t) 

x 2 (t) = G(p)x 2 (t) * ✓Qn(w i )cos( a , 1 t+e 1 ) + n d2 (t) (41) 


where 


- | G ($r)| 


H[ Wi /2ir)/2B s /ii)]' 


(42) 


Bg is the one-sided noise bandwidth of G, Q is the maximum power out 
of the TOR satellite, u»^/2ir and are chosen with uniform pdf's on 
[-10 fHz, 10 MHz] and [0,2*], respectively, and n^ and n d2 are the 
quadrature conponents of the downlink noise. 10 MHz is half the 
bandwidth of the TORS input filter. 

The contribution of the CW pulses to S'(0) is given by 


[S' (0) Jin a pulse of CW RFI] = 0 


(43) 


This is because during such a pulse of a tone characterized by 
sod 0^ , 




-Qn(w-) 


sin toiTrv 

cos(w-T n +2e. ) 

“i T 0 1 0 1 


(44) 


where Tq - 5 usee is pulse duration. Averaging over «. and yields 
zero. 


Now let us obtain Ng conditioned on the presence of CW RFI. 

Nq is the one-sided psd of n, which is the same as the z process 
when = 0. During the ith pulse n is given by 

n(t) = -XjftJx^t) = -Qn(w^) ^ sin 2(u.t+e i ) (45) 

where we have neglected downlink noise. Since Tq « 1/2B^ (the 
correlation time of the H filter), then n during the pulse looks like 
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a 6 -function Input to H of area -Qn(<u.) -Sf- 


tq sin Mii, 


i 0 


The first zero crossing of 


sin w*tq 

~Vo 


“i T 0 


sin(w.To+26i). 

is at u-/2w = 1/2 tq = .1 MHz, 

while n(w^) is still close to one for u^/Iir several times greater 

than that. So the operation of integrating over a pulse duration 

passes only a small fraction of the tones that G does. If we assume 

that the H filter integrates over intervals of length 1/2B L and 

multiplies the result by 2B, , then the output of H due to CW RFI 
a/2B. l ~*- 


x n a/ to i 

is -2 B l Q-f g L n^) 


S1n “1*0 
w i T 0 


s1n(<d. Tq+ 20^ ) where A is the 
repetition rate of the CW RFI process. Since A » we can assume 
that in any interval of length 1/2B L the distribution of the u>^ *s and 
e^'s of the tones is about the same as the statistical distribution 
of and 8j, say. So the variance of that part of H(p)n(t) due 
to CW RFI is given by 


Var 


* 

O T ft F Sin w,T n T 

= B l At 0 Q [ n(w l J -S^ sintY^JJ 


where wj/2it is uniformly distributed on [-10 Miz, 10 MHz], But we 
have -x 

i2 


r f / \ $in “l T 0l . 1 1 

E Ul U)-J Tq ' 2 * 2t ( 


1 


10 MHz 


(47) 


Therefore, the conditioned Nq is 


[Nq | in pulse of CW RFI] = Q 2 . * 


1 


10 MHz 


(48) 


In the evaluations of phase error for S-band done in Section 4.6.3, 
the contribution to the sum Nq made by CW RFI is an order of magnitude 
greater than that made by the more powerful of the two WGN RFI 
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processes and on the same order as that made by the no-RFI situation, 
which has much larger duty cycle. 

4.5.4 RFI Effect on Symbol Tracking 

4. 5.4.1 General De scription of RFI in the loop 

We need to obtain a model for the RFI in the digital -transition 
tracking loop (DTTL) shown in Figure 3.5 and then a new expression 
for timing-error variance. The loop handles the stream of biphase 
symbols as if it were a stream of double-rate NRZ symbols. Just as 
for the carrier loop, the modeling for RFI effects depends on the 
particular S-band symbol rates and RFI environment. The relationships 
among the important parameters aren't as favorable here as for the 
carrier loop, so this model is probably not as good as that one. The 
parameters we consider at this point are two: T, the duration of 
half a biphase symbol which has value 1.736 ysec (mode 1) or .868 
ysec (mode 2); and pulse duration, which is 5 ysec for almost all the 
more powerful RFI pulses. 

In order to analyze this loop in a way similar to the carrier 
loop, we make this loop be continuous-time by considering step 
functions I, J, and e which take on, respectively, the values I fc , J fc , 
and e^ (shown in Figure 3.5) on intervals of length T. Then e is just 
the dynamic phase error process, like the z process for the Costas loop. 
The correlation times of I and J are, respectively, 2T and T, so that 
the correlation time of e lies between those values. Now, 2T = 3.5 
ysec or 1.7 ysec. We will assume, not with great accuracy, that 
pulse durations are enough greater than the correlation time of e 
so that e has two types of characteristics; one during no RFI 
or a pulse of WGN RFI and the other during a pulse of CW RFI. 
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Then (40) holds for S'(0) and so does a similar equation for 
Nq. S * (0) = g'(0) and Nq « 2S(0,0) in the notation of Section 3.5. 

We must obtain the conditioned quantities. 

4. 5.4.2 S‘(0) and Nq During No RFI or WGN RFI 

Just as for the carrier loop, the situation during a WGN-RFI 
pulse is like that of no RFI except for an increase in uplink thermal 
noise. So the conditioned values of S'(0) and of Nq are just obtained 
from (3-76) and (3-79) using conditioned values of symbol r.m.s. 
voltage A and one-sided noise psd Nq. It should be noted that "T" 
in those equations refers to biphase symbol duration. 

4. 5.4. 3 S‘(0) and Nq During CM RFI 

Just as for the carrier loop, we combine the CW RFI processes 
into one with pulse duration = 5 psec and say that during the 
ith pulse the loop input is merely a tone *4) cos(w.t+e.. ), where Q 
is the maximum power out of the TDR satellite and o^/2tt and are 
chosen with uniform pdf's over [-10 MHz, 10 MHz] and [0,2ir], respectively. 
Then the important loop quantities have the values given below: 

1 r / sin(o).T/2) v*| , 

‘k * 2 [ Si9n (-^T72 Jplgn cos( “i (k ‘ jJT+0,) 

- sign cos(k+ ^-IT+e . )3 (49) 

sin(C n w i T/2) 

J k ■ h* W c “("i kT+e i) ( 50 > 

e k * >k J k < 5, > 

The contribution to g'(0) is zero for the CW RFI. 

Ue must find an equivalent Nq conditioned on the occurrence 
of CW RFI. It must be recalled that in the linearized CW phase- 
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locked loop equivalent to the DTTL, a filter H filters the scaled 

equivalent noise process n to yield timing error. The n process 

equals the e process when timing error is zero. Since « 1/2B L , 

where B. is the one-sided loop bandwidth (of H), then the ith 
L M 

pulse looks like a 6 -function input to H of area T ]£ e. , where 

k=l k 

M = t/T = 2.9 (mode 1) or 5.8 (mode 2). If H corresponds to 

integration over an interval of length 1/2B L and multiplication by 

2B, , then the variance of the filter output due to the CW RFI is 

2B, AT E a ( Y, e. ) , where A is the CW pulse repetition rate, 

11 k=l 2 t z , M \ z 

so that the equivalent N« value is E . ( 53 e b J . 

0 T 0 “i * e i ' k^l k/ 

In the evaluations of r.m.s. timing error for the S-band done 
in Section 4.6.3, the contribution to the sum made by CW RFI 
is about the same as that made by the ’re powerful of the two 
WGN RFI processes. 

4 . 6 Preliminary LinCsim Performance Predictions 
4.6.1 Introduction 


Sections 4.6.1 and 4.6.2 have been omitted pending 
clarification of NASA Headquarters TWX 1/0083 of 2 January 
1980 which pertained to the handling of information related 
to TDRSS vulnerability information analysis, data, etc. 
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4.6.3 Predicted Synchronization Performance 


Parts of this section have been omitted pending 
clarification of NASA Headquarters TWX 1/0083 of 
2 January 1980 which pertained to the handling of 
information related to kDRSS vulnerability information 
analysis, data, etc. 


This section gives the predicted . performance of the ground-station carrier/ 
subcarrier phase and symbol -timing recovery for the Shuttle Ku-band 
mode 1 and S-band return links. Results were obtained from analytical 
models implemented as computer programs for the two-channel Costas 
loop with hard-limiters, single-channel Costas loop, and data transition 
tracking loop. Plots of the r.m.s. phase error due to noise are given 
for the carrier and square-wave subcarrier of Ku-band mode 1 and for. 
the S-band carrier for the allowable range of data rates and symbol 
types (NRZ or biphase). Plots of the r.m.s. timing error due to noise 
are given for all three channels of Ku-band mode 1 and the one S-band 
channel for the allowable range of data rates and symbol types. For 
S-band, results are shown both with and without the RFI model given 
in Table4.1. Tables 4. 5,4.6 and 4.7 list the link parameter values.assumed. 

One conclusion that can be drawn from the plots is that for Ku- 
band subcarrier recovery, NRZ and biphase symbols in channel 2 give 
.practically the same results, for any data rate. 
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Table 4i5 .LinCsim Parameter Values for Figures 4.16 through 420 . 

Link: Ku-band mode 1 dedicated return link 
Signal Modulation: UQPSK 
Power Split (I/Q): 4/1 
I Channel (Representing Channel 3): 

Data rate = 2 Mbps to 50 Kbps, NRZ symbols, 
unspread, convolutionally encoded at rate 1/2 
Q Channel: Modulated square-wave subcarrier 
TDRS Front-End Filter Bandwidth: 225 MHz 
TDRS TWT : max AM/PM = 10° /dB, input backoff=10 dB 
Eb/Ho on Downlink: _> 30 dB 
Receiver Carrier-Recovery Loop: 

Two-channel loop with hard limiters and amplitude 
ratio 2/1; arm filters noise bandwidth = § x x 
max{ 2x1 -channel data rate, 22.666 MHz); one-sided loop 
bandwidth = 1 KHz 

Receiver Symbol-Timing Recovery Loop: 

Data transition tracking loop with = .25 and with 

(loop bandwidth)/(2xl-channel data rate) 

.267%, 2 Mbps £ I-channel data rate <_ 37.5 Mbps 

200 KHz/(2xI-channel data rate), 

37.5 Mbps < I-channel data rate 



- 172 - 








fi.M.S. Phne Error (Dog) 



- 175 - 




R.H.S. Timing Error (Symbol Fraction) 











Table 4.6 . Li nCsim Parameter Values for Figures 4 £i through 4.32. 


Link: Ku-band mode 1 dedicated return link 
Signal Modulation: UQPSK 
Power Split (I/Q): 4/1 
I Channel (Representing Channel 2): 

Data rate = 16 Kbps to 2 Mbps (NRZ) or 16 Kbps to 1024 Kbps 
(biphase), NRZ or biphase symbols, unspread, uncoded 
Q Channel (Representing Channel 1): 

Data rate =192 Kbps, biphase symbols, unspread, uncoded 
TDRS Front-End Filter Bandwidth: 225 MHz 
TDRS TWT: max AM/PM = 10°/dB, input backoff = 10 dB 
E b /N 0 on Downlink: _> 30 dB 

Signal Power is Reduced in Receiver (by multiplying it by 8/ w ^) 

to Reflect Effect of Demodulating the Square-Wave Subcarrier 

with a Sine Wave. 

Receiver Carrier-Recovery Loop: 

Two-channel loop with hard-limiters and amplitude ratio 2/1; 

.3, 300 Kbps < R < 600 Kbps 
— max — 

arm-filters noise b.w. = \ x .6, 600 Kbps < R <1.5 Mbps 
in MHz 2 K - max - 

1.2, 1.5 Mbps <_ R max <_ 3 Mbps 

| max(384 Kbps, I-channel data rate}, I NRZ 
max{384 Kbps, 2xl-channel data rate}, I biphase; 
one-sided loop b.w. = 1 KHz 
Receiver Symbol -Timing Recovery Locp: 

Data transition tracking loop with = .25 and with loop 
bandwidth/data rate = .01% for channels I and Q. 
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R.M.S. Tinting Error (Symbol Traction) 
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APPENDIX 

A DERIVATION U EQ. (21) 

Here we derive the eq. (21) of the t» j xt. Suppose that we observe 
the stream of received data over T sec during which period there are 
N * TR symbols. The average error probability over N symbols can 


be written as 


p.(M * r 


e' s 




+ 


(A.l) 


(A. 2) 


where 

W - «s - E - E 

c g 1 

Note that we have assumed Nq(N $ ) > 0 and also allowed to have non- 
integer values for the ratios (t. /T.). In practice, it is always 

K fc a 

possible to choose T sufficiently large so that Mq(N s ) is larger than 

0. The noninteger (i. /T.) implies that the symbol partially hit by 

a RFI pulse is counted as the symbol error caused by the RFI pulse 

weighted with the fractional part of the ratio (t. /T,). This will 

K «. 0 

yield a slightly pessimistic result when th-' ratio (x. /T.) is near 
unity. And it becomes insignificant for (x^ /T d ) > 1. Rewriting 
eq. (A.l) we have 

K c /"k\ 

P e (N s> * 1 + £ l-r) T k ' p e< k c> - P e<°» 

k c M T 
c N. .n , 1 


2 OK 


V 


(P e (kj) - V 0 ’) 


(A.3) 
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By the assumption that the variables n£ and n£ are random processes 

c g 

with a Poisson distribution 


<X,T) n 

P|(n) = — ui — expt-XjT) 


(A.4) 


where 


X^ * pulse repetition rate of i— group of the RFI 
T = observation period. 

He average the p e (N $ ) over the statistics of n£ and nj[ 


W - 1 + E Vk < P e< k c>- P e<°>> 


k c =l 


c c 
n 9 

+ E 

V 


(A.S) 


which is independent of T. Implicitly, we have assumed that T is 

sufficiently large such that the statistics of the variables n^ 

T ^ 

and n^ satisfy the eq. (A.4). Therefore, we conclude that 


lim W 

T-**> 

N 


N N 

E p k ’k P e( k c> + E P k \W * p o P e«°> < A ' 6 > 

: =1 c c k *1 9 9 


wnere 


N 


p 0 * 1 ■ E p k'k ' E P kjk 

k =1 c c k =1 c c 

c c 

and also that eqs. (22) and (23) of the text hold. 


(A.7) 
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5. PROPOSED ESTL RFI TESTING 


5.1 Introduction 

Since NASA/ JSC has a complete hardware simulation system for 
the non-RFI Shuttle/TDRSS Link it is certainly worthwhile to 
extend this simulation capability to RFI environments by adding 
an RFI test generator. For the design of such an RFI test generator 
several things must be taken into account. First the available data 
on the RFI sources have to be reviewed and all useful statistical data 
must be extracted. Then the link performance sensitivity to these 
RFI features must be estimated in order to decide which statistical 
parameters must be modeled in the simulator for reasonably accurate 
performance predictions. The third consideration would be the 
complexity of the resulting simulator. 

. Section 5.2 reviews some of the features which might be found in 
the RFI signal and which, from our experience with LinCsim may 
substantially affect the overall link performance. Section 5.3 presents 
a simple RFI test generator implementation based solely on the statistical 
information contained in the RFI test environments. It is in our 
opinion the simplest test generator which might still give meaningful 
results. 

5.2 Desirable Features for RFI Simulator 

The RFI test generator should duplicate all those known RFI 
characteristics which affect the overall Shuttle/TDRSS link perform- 
ance. This includes the following features: 

(1) Type of interference. There is a significant difference in 

the effect of wideband Gaussian noise and pulsed Q*/ interference 
on BER performance, even if both have the same power. Therefore, 
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both interference types must be generated. 

(2) Power level. The wide range of power levels could be approximated 
by the histogram defined by the unclassified RFI test environments 
discussed in Chapter 4. The realism of the simulation can be 
improved however by approximating the actual distributions with 
steps smaller than 10 dB. 

(3) Pulse duration. The actual pulse length distribution can be 
approximated b, the two values used in the RFI test environment or, 
preferably, by a discrete distribution over more sample points. 

(4) CW frequency. The p.d.f. of the center frequency of the CW 
pulses should be known and modeled. 

(5) Periodicity. If there is any correlation between the RFI pulse 
arrival times it should be modeled if the correlated •'ulses are 
separated by less than the demodulator memory time. 

(6) (Modulation. If the radar pulses have signatures such as chirp 

or phase coding their modeling as CW pulses may be very pessimistic. 
In certain cases it might be more appropriate to use a narrowband 
Gaussian noise to simulate such modulated pulses. 

(7) Pulse overlaps. The simultaneous presence of more than one CW 
pulse might t>e worth simulating because of the intermodulation 
effects. However, the probability of such overlaps should be 
compared « accuracy of the pulse rate estimates before such 
a decision is made. 

5 . 3 Proposed ESTL RFI Simulator 

The RFI simulator design to be proposed is based on the consider- 
ations listed in the previous section, on the RFI test environment 

discussed in Section 4.2 and on the assumptions listed there. 
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For the RFI generator design it was decided that no effort 
should be made to model the overlap of pulses. This is based on 
the following reasons: If pulses with different power levels interact, 
the resulting intermodulation levels are so much smaller than the 
larger of the two signals that they may be neglected. The overlap of 
two pulses of the same power level on the other hand has negligible 
probability. 

The proposed RFI test set up is shown in Fig. 5.1. The RFI test 
generator produces an S-band output signal whose bandwidth and center 
frequency agree with the values of the Shuttle S-band link under test. 
This signal is added to the return link signal between the spaceloss 
simulator the " r o simulator. 

The R I lust ,m<. >.«.or shown in more detail in Fig. 5.2. It 
contains throe cheru^iS, one to generate noise bursts and two to 
generate CW bursts. The gating and power level selection is done by 
a programmable attenuator under control of the burst timing and power 
control logic. 

The frequency of the CW tone bursts is given by the instantaneous 
frequency of the sweep generators. To avoid generating chirp signals 
the sweep rate has to be slow enough to leave the frequency approximately 
constant over one burst. This means however, that consecutive CW 
pulses have very similar center frequencies. The effect of this 
unwanted correlation is minimized by using two sweep generators, one 
covering the data bandwidth, the other the rest of the channel bandwidth. 
Since the inband sweep generator is gated at a very low duty cycle 
it can change its frequency by an adequate amount between pulses to 
minimize the unwanted correlation. 









'cjCinC^Oi 


NOISE 

GENERATOR^| 


BURST TIMING 
+ POWER CONTROL 


| PROGRAMMABLE 
ATTENUATOR j 



Figure 5.2. Proposed RFI Test Generator. 


'dlinCoi 




inL^om 


6. MISCELLANEOUS ANALYSES AMD RECOWgNOATIQNS 

6.1 Introduction 

This chapter documents several analyses performed under the 
current contract. These studies were short-term efforts in support 
of JSC but not directly related to the contract task statements. 

The findings and recommendations were communicated to JSC personnel 
in the form of short memos. 

6.2 Waveguide Effects on PSK Signals 

6.2.1 Introduction 

In a typical satellite terminal there maybealong waveguide 
section from the low-noise amplifier near the antenna to the 
remaining signal processing equipment. Based on measurements taken 
at KSC it was predicted that the linear distortion in a 100 m 
section of WR62 copper waveguide could result in more than 10 dB 
performance degradation. LinCom predicted the degradation analytically 
using LinCsim and found it to be a fraction of a dB. A similar 
analysis was then made for the ESTL/SAIL test setup and an interpreta- 
tion of the KSC measurements was found. 

6.2.2 Linear Distortion of a Waveguide 

The major linear distortion in a long waveguide section is the 
quadratic phase nonlinearity. Fig. 6.1 shows the phase characteristic 
of a 100 m section of waveguide WR62 for the center frequency 15.0034 
GHz. The resulting performance degradation for a 50 Mb/sec BPSK NRZ 
bitstream is shown in Fig. 6.2. 

The ESTL Shuttle/TDRSS link simulator uses a 275 ft section of 
the same waveguide and its effect was also analyzed using LinCsim. 

The results depend on the center frequency of the signal. For the 
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Figure 6.3. Effect of Waveguide on BER at 15 GHz. 
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Figure 6.4. Effect of Waveguide on BER at 14 GHz. 


0037 





SSO-to-TDRS link (f c B 15.0085 GHz) the degradation Is approximately 

.1 dB as shown in Fig. 6.3. For the TDRS-to-Ground link (f * 

c 

13.9375 GHz) the CNR loss is close to .4 dB (Fig. 6.4). This is 
due to the increased linear distortion in closer proximity to the 
cutoff -frequency fg = 9.49 GHz. 

6.2.3 Interpretation of KSC Measurements 

The system assumed for the above computations is sketched in 
Fig. 6.5. There are two loss components which will be observed in 
the measurement setup, viz., signal attenuation and phase dispersion. 
Note that the signal -to-noise ratio is set by the LNA before the 
waveguide run. The signal distortion effect which the waveguide 
produced on E^/Ng is a linear phase dispersion which degrades link 
performance by .2 dB for 100 m waveguide in a coded system and may 
be .5 dB in an uncoded system. If the LNA at the antenna output is 
not present the signal -to-noise ratio is set after the waveguide 
and the waveguide attenuation shows up as an E^/Nq loss. This loss 
is tabulated in Table 6.1. From the magnitude of this loss it is 
obvious that no system should be implemented in this fashion. Also, 
the problem could not be solved by using a different type of waveguide. 

The two sources of loss have different effects on the shape of the 
BER vs CNR curve and it is recommended that the measurement results 
of KSC be analyzed in order to find the source of the losses. If 
the BER curve looks like a shifted ideal BPSK curve of Fig. 6.6, 
the source of the loss is pure attenuation ! If the curve levels 
off as illustrated in Fig. 6.7, the source of loss is phase dispersion . 

6.3 S-Band Low Power Mode 
6.3.1 Introduction 

There is a requirement for the Shuttle S-band link to operate with 
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Figure 6.5. Receiver Configuration Assumed. 



Figure 6j6.B£R Curve for Attenuation Loss. 
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Table 6.1. Waveguide Attenuation. 
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reduced prime power consumption during periods not requiring the trans- 
mission of high data rates. This section addresses the options available 
to reach this goal with minor modifications of the present hardware. 

First the options and the functional requirements of the relevant TDRSS 
services are discussed, then the available Shuttle equipment is reviewed 
and the necessary modifications are identified. Finally a solution is 
recommended . 

6.3.2 TDRSS Services and Functional Requirements 

For S-band return links the TDRSS offers the SSA and MA services. 
However, MA requires a greater user EIRP for a given data rate without 
providing any compensating advantages, hence its use is not recommended 
if the TDRSS SSA service is available. Tables 6.2 and 6.3 list the 
pertinent SSA return link characteristics for PN-spread and unspread 
signals, respectively. It must be noted that during acquisition the 
data modulation must be inhibited in both cases. Also, despite the fact 
that the acquisition EIRP is presently specified 10 dB higher for an 
unspread signal than for a PN-spread one, there is no reason to 
expect the hardware to show this difference in performance. Hence 
the unspread signal can be expected to be acquired with the same 
user EIRP as the PN-spread signal. 

6.3.3 Shuttle Equipment Options 

The Shuttle uses for its direct links to GSTDN, in addition to 
the regular SSO-TDRSS configuration, a low power mode of the above 
system and an FM link, both equipped with an omni antenna. 

6. 3. 3.1 Low Power PM Mode 

Making the low power PM mode (output power 2W, EIRP 1 dBW) 
available for SSO-TDRSS links would require PN spreading to meet the 
acquisition EIRP requirements. However, since there is no justification 
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Table 6.2. SSA Return Link Characteristics for PN-Spread Signals 


Table 6. 


Minimum Data Rate 100 b/s 
EIRP for Above Rate -15 dBW 
Min. Acquisition EIRP -6 dBW 


. SSA Return Link Characteristics for Unspread Signals 

Minimum Data Rate 1 kb/s 
EIRP for Above Rate -6 dBW 
Min. Acquisition EIRP 4 dBW 



for the restrictive specification on the unspread signal it should be 
possible to have the required acquisition EIRP specification reduced to 
less than 1 dBW. Otherwise, the PN modulation would have to be included. 
Since the design 1-dB bandwidth is greater then 5 MHz, it should be 
sufficient for the spread signal. 

Since the resulting signal satisfies regular SSA specifications, 
except for the TEC code used, the impact on the ground station is 
minimal: the soft-decision symbols from the low-rate demodulator have 
to be routed to the rate 1/3 Viterbi decoder. 

6. 3. 3. 2 FM Equi pment 

The FM link (output power 10W, EIRP 4.1 dBW) satisfies the present 
TDRSS EIRP specification for signals without PN spreading. However, a 
phase-modulator would have to be provided in order to obtain a signaling 
format compatible with the ground station equipment. Otherwise, the 
TDRSS IF service can be used with FM demodulation performed by NASA. 

This latter approach may suffer from the fact that the IF channel 
characteristics are presently not well defined due to the RFI 
hardware changes. 

6.3.4 Recommendations 

In the light of the above discussion we propose to use the low 
power PM mode for low rate data transmission. A discussion with Goddard 
personnel about the minimum acquisition EIRP for unspread links is 
recommended before the decision is made to add PN spreading to this 

link. 

6 . 4 Effect of Spacelab Da ta Transition Density on Clock Recovery an d BER 
6.4.1 Introduction 

The Spacelab data do not meet the specifications imposed on the 
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transition density and on the maximum length of transition-free runs. 
The effect of these two factors on the quality of the recovered clock 
and on the bit error probability was estimated using a worst-case 
analysis. The results show that for a reasonable loop bandwidth 
(.IX of symbol rate) the tracking loop can track well through the 
transition free runs and that the recovered clock and the bit error 
rate are only slightly degraded by the low transition density. The 
error rate degradation is expected to be less than .1 dB. The follow- 
ing sections summarize the known facts about the data transition 
statistics and the analysis supporting the above results. 

6.4.2 Characterization of Spacelab Data Transition Density 

Maximum all-zero or all-one runs into t.icoder are 128 bits long. 

Maximum runs occur frequently; up to one every 256 bits. 

Average density is 300 transitions in 1536 bits. 

After the convolutional encoder the maximum all -zero runs will be 
approximately 256 syidbols, spaced 256 symbols apart. The average 
transition density will be more than (300/1536), most likely in the 
range 30* to 40S. 

6.4.3 Analysi s 

The above transition statistics can affect the clock recovery 
loop performance through three mechanisms: 1) the recovered clock 

frequen~y may drift off during transition-free periods due to loop 
stress, 2) the noise-induced phase jitter may increase during such 
periods since the restoring force is missing, 3) the loop signal-to- 
noise ratio is lowered if the transition density is less than 50*. 

If a second-order tracking system is used the first effect may 
be neglected since the drift time-constant is the time-constant of 
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the (imperfect) integrator which is typically many times higher 
than the inverse of the bandwidth. 

The second effect may be upper-bounded by assuming a first-order 
tracking loop. Then the phase error is described by a first order 
Gauss-Markov process during the transition-free period and the variance 
increases by a factor of 4 over 51c symbol -times (neglecting the 
reduction of the variance which takes place between two 256-symbol 
all -zero strings). 

The effect of the low transition density is to reduce the 
effective loop SNR by approximately 2 dB. Since the r.m.s. tracking 
error for a transition density of 50% is epxected to be 1.4% the 
resulting tracking error is 1.7%, yielding in an incremental CNR 
loss of .92 dB [23- Adding the effect of long all zero strings results 
in 3.5% r.m.s. tracking error or approximately .1 dB CNR loss for the 
bit error rate. 

6.5 Spacelab Risetime and Jitter Specifications 

6.5.1 Introduction 

The predicted Spacelab symbol rise-time was increased from 3.5 ns 
to 4.5 ns and the combined bit jitter and data asymmetry from 1 ns to 
2 ns. In this secti >n the effect of these changes on the overall link 
performance is computed. 

6.5.2 Analysis 

The basic pulse shape for the highest data rate is sketched in 
Fig. 6.8. Due to data asymmetry and bit jitter the actual transitions 
are within +1 ns from their expected time (Fig. 6.9). Since the 
scrubber implementation assures that the sampling point is at least 
5 ns from the worst case transition [1] it occurs always beyond the 
90%-point of the transition and the signal-to-noise ratio is very 
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little degraded. The remaining 8 ns interval allows approximately 
6% r.m.s. jitter on the recovered clock while the predicted value is 
less than 22 [2]. This value increases only slightly for increased 
bit jitter since it is dominated by thermal effects [2,p.l26]. He 
may therefore conclude that the overall link performance is not 


measurably degraded. 
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